Structure of symmetric and asymmetric "ripple" phases in lipid bilayers 
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We reproduce the symmetric and asymmetric "rippled" Ppi states of lipid membranes by Monte 
Carlo simulations of a coarse-grained molecular model for lipid-solvent mixtures. The structure and 
properties compare favorably with experiments. The asymmetric ripple state is characterized by a 
periodic array of fully interdigitated "defect" lines. The symmetric ripple state maintains a bilayer 
structure, but is otherwise structurally similar. The main force driving the formation of both ripple 
states is the propensity of lipid molecules with large head groups to exhibit splay. 

PACS numbers: 87.16.Dg, 87.16.Ac, 82.70.Uv, 87.14. Co 



Membranes are ubiquitous in all living organisms 
Their central structural element is a lipid bilayer, which is 
stabilized by the amphiphilic character of lipid molecules 
- they self-assemble such that their hydrophilic head 
groups shield the hydrophobic tails from the surround- 
ing water. Pure lipid bilayers have been studied for 
a long time as model systems that can provide insight 
into the structural properties of biomembranes. Already 
these seemingly simple systems exhibit a rich spectrum 
of structures and phase transitions [1, 0, 0, [H| . The most 
common state in nature is the so-called "fluid" state 
which is characterized by a large number of chain defects 
and high lipid mobility. If one decreases the temperature, 
one encounters a phase transition (the "main" transition) 
to a "gel" state where the lipid molecules are more ordered 
and less mobile. The structure of this low temperature 
phase depends on the interactions between the lipid head 
groups. Loosely speaking, lipids with small head groups 
such as phosphatidylethanolamines Q assume a state 
where the long axes of the chains remain perpendicular 
to the bilayer normal {Lp phase). Lipids with large head 
groups and relatively strong head-head attraction such as 
phosphatidylcholines [H| exhibit tilt {Lp' phase). Finally, 
lipids with large head groups and weak head-head at- 
traction such as as ether linked phospatidylcholines 0, 0| 
form a phase L"p where both opposing lipid layers are 
fully interdigitated 

The main transition has attracted considerable inter- 
est, since it occurs at temperatures that are typical on 
earth (between — 20°C and 60°C). The mechanism that 
governs the transition La Lp to the untilted gel is 
comparatively straightforward. The transition is driven 
by the competition of the entropy of chain disorder and 
the free energy of chain alignment 0, @| (i.e., chain pack- 
ing) and is thus in some sense related to the isotropic- 
nematic transition of liquid crystals. At the transition 
La ^ Lfji to the tilted gel, the situation is much more 
complicated. Here, the main transition is preempted by 
a "pretransition", and one observes an intermediate state 
with a periodic, wave-like surface structure: The "rip- 
ple" phase Pfj' , first reported by Tardieu et al. The 
microscopic structure of this mysterious phase has been 



debated for a long time. 

In fact, at least two different rippled states have been 
reported, which often coexist Electron density maps 
(EDMs) have recently been derived for both of them from 



X-ray scattering data [Ul, |l^ . One of the structures is 



asymmetric and has a sawtooth profile with alternating 
thin and thick arms and a periodicity of 13-15 nms, which 
corresponds to roughly 20 lipid diameters. The other one 
is symmetric and has a wavy sinusoidal structure with 
twice the period of the asymmetric structure [l3|. The 
formation of the ripples depends strongly on the thermal 
history 13, 13, 3- If the membrane is heated up from 
the gel state, asymmetric ripples are obtained. If one 
cools down from the fluid state, both types of ripples are 
formed - predominantly asymmetric ones if the cooling 
was fast, and predominantly symmetric ones if the cool- 
ing was slow and if a long time was spent at the transition 
temperature. Dynamical X-ray [3| and AFM [l3| stud- 
ies suggest that the symmetric ripple state is metastable 
and very slowly transforms into the asymmetric ripple 
state; however, this does not yet seem to be fully settled. 
The degree of ordering in the ripple states largely resem- 
bles that in the gel state, except for a certain amount of 
disorder in the structure j3| - calorimetric studies sug- 
gest that approximately 10 % of all chains are melted. 
Most strikingly, the self-diffusion of lipids in the ripple 
states is a few orders of magnitude higher than that in 
the gel state, and highly anisotropic [1^|. This has lead 
to the assumption that the ripple states might contain 
"coexisting" gel-state and fluid-state lipids. 

Numerous theoretical models for the ripple phase have 
been proposed, which explain the ripple formation by dif- 
ferent mechanisms: Chain packing effects 17. 18|. dipolar 
interactions [l^l, a coup lin g of monolayer curvature with 
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22 1 or with tilt 
24| . This list is far from com- 



local chain melting (IE [III, 

combination with chirality 
plete. In contrast, molecular simulations of rippled mem- 
brane states are still scarce. Kranenburg et al. (i^l were 
the flrst to reproduce a periodically modulated mem- 
brane state in a dissipative-particle dynamics (DPD) sim- 
ulation of a coarse-grained lipid model. They observe a 
periodic sequence of stripes with alternating gel and liq- 



2 



uid order, similar to a structure proposed theoretically 
by Falkovitz et al. Unfortunately, the distribution 

of head groups in that structure is not consistent with the 
experimental EDMs - the structure is neither asymmet- 
ric, nor does it feature the waviness which characterizes 
the symmetric ripple. Moreover, the relative fraction of 
molten molecules - 50 % - seems too high, compared to 
experiments. A second, very interesting simulation has 
recently been carried out by de Vries et al. jH]. In an 
atomistic model of a lecithin bilayer, these authors found 
a structure containing a stretch of interdigitated mem- 
brane and a stretch of gel membrane. The interdigitated 
patch connects the neighboring gel membrane such that 
the upper leaflet of the bilayer on one side crosses over 
into the lower leaflet on the other side. The authors as- 
sume that this structure will repeat itself periodically in 
larger systems and identify it with an asymmetric ripple. 
It is worth noting that the lipids are not arranged in a 
continuous bilayer - as had been assumed in all previous 
models for the ripple state. 

In this letter, we present Monte Carlo simulations of 
a simplified coarse-grained lipid model, which reproduce 
asymmetric and symmetric ripple states with properties 
that compare very favorably to experiments. The struc- 
ture of the asymmetric ripple is similar to that proposed 
by de Vries et al. [i^. Our simulations show that it is 
indeed a periodic structure, and that it is generic, i.e., it 
does not depend on molecular details of the Hpids. More- 
over, they enable us to propose a structural model for the 
symmetric ripple as well, and to identify the mechanisms 
that stabilize the rippled structures. 

We employ a Hpid model which we have used earlier 
to investigate phase transitions in Langmuir monolay- 
ers [irl . [iil : Lipid molecules are represented by chains 
made of one head bead and six tail beads (Fig.[T}3), which 
are connected by anharmonic springs and subject to an 
intramolecular bending potential. The tail beads attract 
one another with a truncated and shifted Lennard-Jones 
potential (diameter a, well depth e). The head beads 
are larger than the tail beads (I.Ict) and purely repulsive. 
The other parameters and the exact form of the poten- 
tials can be found in Ref. [H (the model corresponding 
to Fig. 7). Self-assembly of the "Hpids" is enforced with 
a recently proposed "phantom solvent" environment [29| : 




FIG. 1: Illustration of our lipid model and snapshots (side- 
view) of two lipid states (at psksT = 2(./a'^). (a) All-atom 
model of DPPC (b) Coarse-grained bead-spring model used 
in this work (c) The fluid phase La at RbT = 1.3e (d) The 
tilted gel phase L^i at fcsT = l.le. For better visualization, 
only heads (reduced size) and tail bonds are shown. 




FIG. 2: Phase diagram of the lipid model as a function of tem- 
perature, T, and phantom solvent density i.e., effective head 
interaction, pa. The phases are: La (fluid), L^/ (tilted gel), 
LJj' (interdigitated gel), (ripple). At high temperatures, 
the bilayer disintegrates. Open squares indicate transition 
points from simulation runs of small systems initially set up 
as fully ordered, untilted bilayers. No ripples were observed in 
that case. Closed squares show transition points determined 
by heating up and cooling down the system (different system 
sizes), with the error bars corresponding to the width of the 
hysteresis. Open circles denote points where the membrane 
undergoes a transition from L^/ to L"' upon heating. 

We add "solvent" particles which interact only with Hpid 
beads (repulsively), and not with one another. Two ex- 
amples of self-assembled membranes in the La and the 
Lp' state are shown in Fig. [T]:,d). The phantom solvent 
has the simple physical interpretation that it probes the 
accessible free volume for solvent particles in the pres- 
ence of lipids. It entropically penalizes lipid/solvent in- 
terfaces, and thus effectively creates an attractive "deple- 
tion" interaction between the lipid beads next to such an 
interface, i.e., the head beads [30||- The strength of this 
interaction is directly proportional to the phantom sol- 
vent density, ps. Compared to other expHcit solvent en- 
vironments, the phantom solvent has the advantage that 
it does not introduce an artificial solvent structure and 
artificial solvent-mediated interactions between periodic 
images of bilayers. Moreover, it is computationally cheap 
- in Monte Carlo simulations, less than 10 % of the to- 
tal computer time is typically spent on the uninteresting 
solvent region. 

We have carried out Monte Carlo simulation at con- 
stant pressure with periodic boundaries in a simulation 
box of fiuctuating size and shape. This ensured that the 
membranes had vanishing surface tension. The system 
sizes ranged from 288 to 1800 lipids (corresponding to 
2000-12600 beads), typical run lengths were 1-10 mil- 
lion Monte Carlo sweeps. The resulting phase diagram is 
shown in Fig. [2l The model reproduces the experimen- 
tally observed gel and fiuid phases for lipids with large 
heads: The fiuid phase {La), the interdigitated gel {L'^*) 
for low Pa, i.e., weak head attraction, and the tilted gel 
(L/3') for higher ps, i.e., strong head attraction. The 
structures of these phases and the phase transitions shall 
be discussed in detail elsewhere [3l|. 

Here, we fix the solvent density at PsksT = 2e/o'^, 
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where the gel phase has the tilted Lpi structure. In the 
transition region between Lpi and La, we observe modu- 
lated configurations which we identify with rippled states. 
They develop spontaneously and reproducibly when cool- 
ing a fluid membrane or heating a tilted gel membrane in 
a temperature range close to the transition temperature. 
As in the experiments, their exact structure depends on 
the thermal history. Fig. [3] (left) shows three examples. 

The structure in Fig.jSK) emerged after cooling the sys- 
tem rapidly from the fluid phase down to ksT = l.le. It 
exhibits two ripples of width ^ Iba with a structure very 
similar to that found by de Vries et al. [i^ in their atom- 
istic simulations: At each ripple, a thin interdigitated line 
defect connects a lower with an upper monolayer. The 
second monolayer ends at that line with an edge of disor- 
dered, melted chains. The period, ^ 15cr, corresponds to 
the "natural" period of asymmetric ripples in our model. 
This was deduced by comparing simulations of systems 
with different sizes. The quality of the ripple formation 
depends on the initial box dimensions (before cooling). 
In systems of half the size, only one ripple formed. If 
one box dimension was close to a multiple of 15tT, the 
box shape readjusted to accomodate the optimal ripple 
width. If the initial box dimensions were very unfavor- 
able (20 or 24 cr), no clean ripples formed; instead, the 
ripples bulged and developed interconnected structures. 

The second structure, shown in Fig. [3b), resulted from 
heating up a bilayer in the tilted gel state up to a temper- 
ature close to the main transition, fc^T = 1.21e. During 
the simulation, the bilayer flrst fluctuated very strongly. 
After ~ 6 miUion Monte Carlo sweeps, the tilt was so 
strong that the lipids in both monolayers sHd along each 
other and connected with the other monolayer. The flnal 
structure exhibited one asymmetric ripple and fluctuated 
much less. This structure remained stable for another 6 
miUion sweeps. Apparently, the formation of the second 
ripple is prevented kinetically. 

The third structure. Fig. [3b), shows a membrane which 
has been cooled down from the fluid state to a temper- 
ature close to the main transition fcsT = 1.18e. In this 
case, a new type of structure emerged: The membrane 
maintains its bilayer structure, but the monolayers con- 
tain curved, ordered stripes with a width of roughly 25 
a. These "gel" stripes on the upper and lower monolay- 
ers are interlocked, such that the membrane assumes an 
overall sinusoidal shape. Each stripe ends on both sides 
with conical regions of disordered chains, which are very 
similar to the monolayer caps in the asymmetric ripple 
state. The total width of a ripple is ^ 30cr, which is 
twice as much as the width of an asymmetric ripple. We 
identify this structure with the symmetric ripple state. 

To support this hypothesis, we superimpose the pro- 
posed structures for the asymmetric and the symmetric 
ripple with the EDMs of Sengupta et al. [l3| in Fig. [H 
They can be inserted very nicely. In the asymmetric case, 
they explain the sawtooth shape with the thin and thick 




FIG. 3: Three examples of ripple configurations in a model 
bilayer of 720 molecules (left) with corresponding tilt profiles 
6(x) (right) (a) Two asymmetric ripples, formed after rapidly 
cooling down to the temperature ksT — l.le from the La 
phase (b) One asymmetric ripple, formed after heating up to 
ksT — 1.21e from the L/ji phase (c) Symmetric ripple, formed 
after slowly cooling down to the temperature UbT — 1.18e 
from the La phase. Note that the snapshots (left) show side- 
views of whole fluctuating conflgurations - therefore, a large 
number of heads seem to be buried inside the layer, even 
though they are at the surface. The curves in the graphs 
(right) are shifted in the x direction and replicated periodi- 
cally. Open circles correspond to the lower monolayer, closed 
squares to the upper monolayer. The hatched ellipses indi- 
cate the regions with the interdigitated line defect, and the 
thick solid lines in c) the slopes of 6 on the ordered monolayer 
regions. 

arm (assuming that the interdigitated region is indeed 
small). In the symmetric case, they reproduce the sinu- 
soidal shape. 

The simulations can be used to characterize the ripple 
states in more detail. We just summarize some of the 
results here, the data will be presented and discussed 
elsewhere [3l|. The structure of the ripple states is in 
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FIG. 4: Sketch of the proposed microscopic structures for the 
ripple states layers superimposed onto EDMs from Ref. [l^ . 

(a) asymmetric ripple (on an EDM for DMPC at 18.2 " C) 

(b) symmetric ripple (on an EDM for DPPC at 39.2 " C) 
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many respect similar to that of a gel: Roughly 85% 
of the chains have chain lengths distributed as in the 
gel, only 15% have a reduced length. This is in rough 
agreement with the experimental findings on the amount 
of chain disorder in the ripple state. The head layer in the 
ordered parts of the ripple state has the same thickness 
than in the gel state. The structure factor of the ripple 
state indicates a large amount of positional order, and 
resembles that of an untilted gel. The most revealing 
structural feature is the average tilt of the molecules. 
It points perpendicular to the ripple and is modulated. 
Fig. [3] (right) shows profiles of the average tilt angle 9 
for the three ripples discussed above. The slope 9'{x) 
turns out to be almost constant throughout the whole 
ordered part of the monolayer. Moreover, the numerical 
values are comparable: 6*' ~ 2.6 /cr on average for the two 
asymmetric ripples, and 6' 2.5/tT for the symmetric 
ripple. Even the single asymmetric ripple of Fig. [SJd), 
which has an unfavorable period, still features a constant 
slope of 6*' r-^ 2.3/ct. This suggests strongly that the ripple 
formation is primarily driven by the propensity of lipids 
with large head groups to exhibit splay. The driving 
force of the local curvature, i.e., head packing, is much 
weaker. It does however play a role in determining the 
periodicity of the ripple. Our cooling simulations clearly 
indicate that the system favors a certain ripple width, for 
which both splay and curvature are optimized. 

In sum, we have reproduced symmetric and asymmet- 
ric rippled states with a generic model for lipid mem- 
branes. The comparison with experiments is favorable: 
The structure is consistent with the available EDMs, the 
period length is of the same order as the experimental pe- 
riod length (~ 15 lipid diameters), the amount of chain 
disorder is comparable, and we observe the same depen- 
dence on thermal history. Therefore, we believe to have 
strong evidence that our structures correspond to the real 
ripple states observed in experiments. Factors that are 
important for the formation of these states are: (i) The 
vicinity to the La phase, such that a small number of 
chains can melt, (ii) a strong tendency of monolayers 
to splay inwards - caused by a mismatch between head 
group and tail size, and (iii) the possibility to interdig- 
itate. Chirality is not necessary, in agreement with ex- 
periments [3| ; the "lipids" do not even have to be asym- 
metric. 

In fact, the factor (iii) is only needed to stabilize the 
asymmetric ripple state. If it is absent, the system can 
still form a symmetric ripple state. We note that the 
symmetric and the asymmetric ripples are structurally 
quite similar. Both contain about the same amount of 
molten chains, both have large ordered monolayer regions 
with comparable splay. This explains why the two types 
of ripples coexist, and why it seems so hard to determine 
which one is stable. Our results suggest that the answer 
to that question may depend on the type of lipid, e.g., 
on the head interactions and other factors that promote 



or prevent interdigitation. Coarse-grained lipid models 
may help to study this systematically. Unfortunately, we 
have not yet been able to develop an efficient strategy 
to determine the free energy difference between the two 
states. This will be subject of future work. 
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EDMs, and the NIC computer center Jiilich for com- 
puter time. This work has been funded by the Deutsche 
Forschungsgemeinschaft within the SFB 613. 
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